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We study the Landau model of the class of incommensurate systems with a scalar order parameter 
where the modulated phase is driven by a gradient-squared term with negative coefficient. For 
example, theoretical studies of cholesteric liquid crystals in a field (electric or magnetic) suggest that 
such an modulated phase should exist at high chirality. The bulk phase diagram in the presence of a 
bulk external field which couples linearly to the order parameter exhibits a modulated phase inside 
a loop in the temperature-field plane, and a homogeneous phase outside. On analyzing the same 
model for a semi-infinite system, we find a surprising result; the system exhibits surface states in a 
region where the bulk phase is homogeneous (but close to the modulated region). These states are 
very different from the well-known surface states induced either by a surface field or by enhanced 
interactions at the surface, for they exist and are energetically favored even when the sole effect of 
the surface is to terminate the bulk, as expressed by free boundary conditions taken at the surface. 
Near the surface, the surface-state order parameter is very different from the bulk value (in fact, it 
has the opposite sign). When the temperature or the bulk field are varied to move away from the 
modulated state, we find a surface phase transition at which the surface states become energetically 
unfavorable, though they continue to exist as metastable states. We then study how a surface field 
changes the surface phase diagram. 



PACS numbers: 61.30.Cz; 64.60. Kw; 64.70.Md 
I. INTRODUCTION 

It is well known that a surface field can give rise to wet- 
ting phenomena and also that enhanced interactions near 
a surface can give rise to surface order without bulk order 
101 . Nakanishi and Fisher |^ have given a unified picture 
of wetting and surface ordering at the phenomenological 
(Landau-theory) level; these effects require that surface 
terms be added to the bulk free energy. In this article, we 
report an entirely new surface effect which should occur 
in a particular class of incommensurate systems. We find 
that surface states exist and are energetically favored by 
the mere presence of the surface, without surface terms 
like those considered in 

Candidate physical systems for observing these states 
include highly chiral cholcstcrics in electric or magnetic 
field, where a bulk undulating phase was recently pre- 
dicted to occur [|. This phase is an undulating structure 
in which the amount of orientational order varies periodi- 
cally in conjunction with an oscillation of the direction of 
the local optic axis. It is expected to occur under appro- 
priate conditions of temperature and a strong aligning 
electric or magnetic field. As discussed in section III C 
of reference ||] , the order parameters for the modulated 
state are the amplitudes of the harmonics. The free en- 
ergy that results is identical to that of Landau models 
in which the coefficient of the gradient-squared elastic 
terms is negative, necessitating the inclusion of terms 
quadratic in second derivatives. When the coefficient of 
the gradient-squared term vanishes, a Lifshitz point oc- 



curs in the phase diagram Q . Therefore other candidates 
include Lifshitz-point systems such as the magnetic ma- 
terial MnP laJq], and Langmuir monolayers and diblock 
copolymers [Qwith modulated phases. 

Our Landau model gives a bulk temperature-field 
phase diagram with a closed loop separating the modu- 
lated phase (favored inside) from the homogeneous phase. 
The surprise is that the mere existence of the surface pro- 
duces a surface state which is energetically favored within 
a second closed loop well outside the first. Outside the 
second loop, the surface state exists but it is metastable 
(the equilibrium solution is simply the homogeneous bulk 
state). The order parameter in the surface state is not 
a small perturbation to the bulk order parameter. The 
width of the state has no pronounced temperature depen- 
dence; in particular, it does not diverge. In the presence 
of a surface field (coupling linearly to the order param- 
eter), the surface phase may still occur but the line of 
surface transitions no longer forms a closed loop. ^From 
the above and other evidence, our surface states are very 
different from the states considered in [1,2]. 

This paper is organized as follows. Section II presents 
the efi^ective Landau-Ginzburg model and then describes 
analytical and numerical results for the bulk phase dia- 
gram. Section III presents analytical and numerical re- 
sults for the surface states, in the absence of a surface 
field. Section IV shows how the surface phase diagram is 
modified by a surface field. Finally section V discusses 
the results and their possible realization. 
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II. MODEL AND BULK PHASE DIAGRAM 
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In this section we introduce the model used in the rest 
of the article, and we study the bulk phase diagram, es- 
pecially the transition line separating the homogeneous 
and modulated states. The bulk free energy F), is the 
spatial integral of the density jFf,, which is the following 
functional of the scalar order parameter 4>{x): 

= -hcb+'^rc^' + - i(<^')' + (1) 

where (j)' = d(f>/dx. We have scaled the order parameter, 
the energy and the unit of length to simplify the coeffi- 
cients, and so h and r arc the rescaled ordering field and 
temperature variables respectively. The corresponding 
Euler-Lagrange equation is 

0"" + <l>" ~ h + r(j> + (l)^ ^0 . (2) 

Nakanishi and Fisher examined a very different model; 
the gradient-squared term appeared with a positive co- 
efficient, the ((/)")^ term was omitted, and surface terms 
were added. Their model, without the surface terms, ap- 
plies to the usual Ising model with ferromagnetic interac- 
tions; it has only the disordered and homogeneous (fer- 
romagnetically ordered) phases, and its bulk (r, h) phase 
diagram consists of a first-order line at ft. = and r < 
which terminates at a critical point at r = 0. The model 
of Equation but without the {(f>')'^ term, exhibits a 
Lifshitz point at /i = r = and a first-order line for 
r < 0. 

Without the bulk field h, the model (|l|) has a dis- 
ordered phase at high temperature (T), a second-order 
transition at r = i to a modulated phase, and a strong 
ffi'st-order transition at r « —1.2 to one of two degener- 
ate homogeneous phases; the modulated phase is almost 
sinusoidal over its entire range, and its wavenumber is 
almost independent of T. In the (r, h) plane, the mod- 
ulated phase occupies a closed loop . Outside this 
loop, the energetically favored phase is the homogeneous 
phase, with order parameter (pQ found from 

- ft + r0o + <^o' = ; (3) 

its free-energy density is J-q — — ft0o + ^'^^o + J^o- 

Figure 1 gives the bulk phase diagram, as found for the 
most part by numerical solution of the Euler-Lagrange 
equation (j^) with periodic boundary conditions. The 
homogeneous-modulated transition is second-order near 
r = i, but otherwise first-order. The second-order seg- 
ment and the tricritical points at its ends are found an- 
alytically in the following. 

We consider a spatially modulated order parameter 
and expand it in harmonics. If q is the wavenumber of the 
modulated structure and e is the amplitude of the leading 
harmonic, then the order parameter takes the form 

(j){x) 00 + e cos{qx) + [<j)2 cos(2ga;) + (^2] + 0{e^) (4) 
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FIG. 1. The {r,h) bulk phase diagram corresponding to 
model (1). The undulating state (U) is energetically preferred 
inside the loop, and the homogeneously ordered state out- 
side it. The transition between the states is either first-order 
(solid line) or second-order (dashed line); two tricitical points 
(solid circles) separate the two types of transition. The 
first-order segment was found by numerical solution of the Eu- 
ler-Lagrange equation, the second-order segment from Equa- 
tion (8), and the tricritical points from Equation (10). At 
lower temperatures, a first-order transition at ft = sepa- 
rates the two homogeneously ordered states 4> > and (f> < 0. 

where 4)2 and (j32 are constants to be determined. In- 
serting this order parameter in the free energy (|l|) and 
integrating over a period, one finds the following expan- 
sion of the free energy (per unit volume): 

(.Ffc) =.Fo + e2^2 + e^^4 + 0(e6) (5) 

with coefficients 

.F4 = i(r + + 2)02 + 1 + 3^2)^2 

+ §0002 + 10002 + ^ . (7) 

The free energy has already been minimized with respect 
to the wavenumber q, giving q = + 0{e^). The 

homogeneous phase is unstable to a modulated pertur- 
bation when T2 < 0. Provided that JF4 > then, a 
second-order transition occurs at 

When .F4 is negative, the transition to the modulated 
phase is first-order. To find the tricritical points separat- 
ing the continuous and first-order segments, we minimize 
JF4 with respect to 02 and 02 and then set the result 
equal to zero. On the line J-2 = 0, JF4 is minimized by 
02 — — 20o/3 and 02 — — 60o; the minimum value is 
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^4 = ^ - f 0^ (9) 
and so the two tricritical points arc located at 

r = 29/152 , h = ±^6/193 . (10) 



III. SURFACE PHASE DIAGRAM 

In this section we consider the surface phase diagram 
of the model for a semi-infinite system, with no sur- 
face field. The presence of the surface generally produces 
states localized near the surface, and the states are ener- 
getically favored in part of the phase diagram. We stud- 
ied the surface states in the region where the bulk phase is 
homogeneous, and examined their transitions with vary- 
ing temperature and the external field. Only a cursory 
examination was made in the region where the bulk is 
modulated; in this region, we found many solutions of 
the Euler-Lagrange equation, so many that a detailed 
analysis was felt unjustified at this time. That is, surface 
states and surface phase transitions may exist inside the 
bulk modulated loop, but have not been studied. 

We consider a system occupying the half-space a; > 0, 
and we assume that the order parameter depends only on 
X. The bulk energy Fi, is found by integrating the density 
of (|lj). In this section, we treat the surface very simply, 
by assuming that it merely terminates the bulk; we thus 
take free boundary conditions at the surface. In section 
IV, however, we assume that the surface also applies a 
local ordering field hg ; then the total energy is Ff, + Fs, 
where 

Fs = ~h,(j), (11) 

and (f)s is the order parameter at a; = 0. The general 
boundary conditions are then 

cj,'^ + C - /i. - , 0^' = . (12) 

We solved the Euler-Lagrange equation ^ numerically 
subject to the boundary conditions ([l^). This equation 
can have many solutions, depending on the bulk field h 
and the temperature variable r. Figure 2 gives the surface 
phase diagram for /ig = 0, as found from examining these 
solutions. The surface states are energetically favorable 
inside the outer loop of the figure (with the qualification 
noted above), and the homogeneous states outside; the 
surface states exist (as solutions of the Euler-Lagrange 
equation) outside this loop but are only metastable there. 
An interesting feature is that the surface orders at r = 1 
for h = 0, but the bulk orders only at r = i. Many more 
surface states were found at lower temperatures, but they 
were always metastable. 

Figure 3 shows a typical profile of the surface state in 
the ordered region, for a small and negative bulk field (to 
break the symmetry) and hs = 0. The order parameter 
decays to the bulk value (which is negative) far from the 
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FIG. 2. The (r, h) surface phase diagram for hs — 0. The 
surface states are energetically favorable inside the outer loop 
(with the qualification noted in the text), and the homoge- 
neous states outside. The transition is first-order except at 
the isolated point r — 1, h — (indicated by a dot) where 
it is continuous. The inner loop is the loop of Figure 1 (the 
scale precludes display of details). The leftmost of the three 
horizontal lines at /i = is the bulk transition between the or- 
dered states (j) > and (p < 0. At the rightmost line (deished), 
the surface state changes discontinuously and the bulk state 
continuously. At the middle line, both surface and bulk states 
change discontinuously. 

surface, but it is large and positive near it; the overshoot- 
ing and the damped oscillations result from a complex 
decay constant, as shown below. Correspondingly, when 
the bulk field is positive, the order parameter of the sur- 
face state is negative near the wall and then decays to the 
positive bulk value. Thus at /i = there is a first-order 
transition at which the surface state changes sign. 

To provide an analytical understanding of these numer- 
ical results and also those of the next section, we present 
the following stability analysis of the homogeneous bulk 
state. The analysis is valid when the deviation of the 
order parameter from the bulk value is small. 

The order parameter is written as (f) = (f)o + 4' , where 
(/)o is given by (^ and Tp is the deviation. The free-energy 
density T = Th ~ associated with i/j is 

T^\{r^ 30o')^' + + - W f + Wf ■ 

(13) 

The energy is minimized by an order parameter -0 which 
satisfies the Euler-Lagrange equation 

(r 4- 300^)^^ + 3000^ + V'^ + V'" + = . (14) 

To prepare for the next section, we include also the sur- 
face free energy ([ll|). The boundary conditions are then 

^-'(0) + 7/^"'(0) - h, = , V"(0) = . (15) 
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FIG. 3. Order parameter 4>{x) of the surface state for pa- 
rameters r = —2, h = 0~ and hs — 0. The order parameter 
is large and positive near the surface; it crosses zero and then 
decays to the bulk value for large x. 



For hs = 0, the homogeneous bulk state = is clearly 
a solution of Equations ( [l^ ) and ([l5|). This solution is 
stable over some region of the (r, h) plane, but it becomes 
unstable at the transition to the bulk modulated state. 

To study the surface states, we solve Equations ([l^ ) 
and ( pi] ) perturbatively in ip. The expansion starts from 
the solution 



Mx)=Ae-°"'+A*e- 



(16) 



of the linearized Equation (|lj). The amplitude A and 
the decay constant a arc both complex; the latter (with 
positive real part) is found from 



i(-l + *7) 



(17) 



where I = and 7 = [4(r-|-3(/)Q) — 1]^/^. The condition 
''/'i(O) = in ([l5| ) gives the amplitude A in terms of 
m — -0(0) as 



A 



m 



1 



(18) 



It is convenient to take m as the expansion parameter. 

The solution ([l6| ) gives the free energy to order m^. 
In order to obtain the free energy to the required order 
(m^), one must find the higher-order contributions to 
-0. Let ip = tpi + 1P2, where 'ip2 is the nonlinear part of 
-0. Inserting this form in Equation (p^, using (16) and 
keeping terms to m^, one finds 

C2e-(2"+"*)^ + ZJe-"^ 4- c.c. (19) 
with ?/'2(0) = ^'2(0) = ^-iid coefficients 
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On using the result -0 = V'l + ^^2 in the free energy (|l3| ) 
and integrating over x, one finds that the free energy of 
the surface state (per unit area) is given by 
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where 



a2m' 
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+ O(m^) (20) 



(21) 



The amplitude of the surface structure is determined by 
minimizing the free energy with respect to m for given 
surface field hg. This amounts to satisfying the first con- 
dition in Equation ( [Tsl ) . We now use the free energy ( |20| ) 
to discuss the surface phase diagram in the region where 
m is small. 

Consider first the case hg = 0. For zero bulk field h, 
00 — and there is no surface state when 02 > (that 
is, i/j = 0). Setting 02 = 0, one finds a continuous tran- 
sition at r = 1 from the disordered bulk state to = 
to a surface state with m ^ 0; this is the second-order 
point at the right of Figure 2. For field h 0, the bulk 
order parameter 0o is also non-zero and the free-energy 
expansion (^ ) has a term in m^; this cubic term gives a 
first-order transition to the surface state, again as found 
numerically. Near the point (r = 1, /i = 0), m is small 
and the transition line can be found approximately from 
the free-energy expansion (20). Away from this point, 
however, the full free energy must be minimized numer- 
ically; Figure 2 gives the resulting (r, h) surface phase 
diagram for hg = 0. 



IV. EFFECT OF A SURFACE FIELD 

We consider now the surface phase diagram for non- 
zero surface field hg. Positive hg, for example, tends to 
increase the order parameters of all states in the region 
near the surface. The new feature is that the Euler- 
Lagrange equation must now be solved numerically for 
what were homogeneous bulk states at ft-s = 0; for lack 
of a better term, we refer to these surface-field-modified 
bulk states simply as bulk states. Figures 4 and 5 give 
parts of typical phase diagrams for hg > 0, as found by 
numerical solution of Equation (|^), subject to the bound- 
ary conditions ([l^). 

Figure 4 shows the high-temperature part of the phase 
diagram for hg = 10^*. The surface field breaks the 
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FIG. 4. The high-temperature part of the (r, h) surface 
phase diagram for hs = 10~*. Both transition hues are 
first-order. The lower line ends at a second-order point 
marked by the solid circle (see text). 

transitions of Figure 2 into two first-order lines at which 
the surface state changes discontinuously. In the region 
bounded by the upper hne and the left vertical (where 
the bulk field h is positive), the order parameter of the 
surface state is negative at the boundary x = {(j)s < 
0). The lower line ends at a second-order point. Below 
this point there is a first-order transition between the 
paramagnetic state and the surface state with cjjs > 0, 
while above it the two states are indistinguishable. The 
free-energy expansion (^^ can be used to find this point 
to leading order in h^] the result is 

in good agreement with the numerical results. 

Figure 5 shows the low-temperature part of the phase 
diagram for hg — 10~^. Paradoxically, a positive surface 
field cooperates, rather than competes, with a negative 
bulk field to enhance the stability of the lower surface 
state (and it competes with a positive bulk field for the 
other). These effects occur because the order parameter 
of the surface state changes sign (as seen in Figure 3). 

V. DISCUSSION 

We have developed and analyzed a model to describe 
the effect of a substrate (or a free surface) on a mate- 
rial which has a bulk phase transition between homo- 
geneous and modulated states. Modulated states tend 
to form because the free energy of the model contains a 
term, quadratic in first derivatives of the order param- 
eter, which has a negative coefficient. We treated the 
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FIG. 5. The low-temperature part of the (r, h) surface 
phase diagram for hs — 10~^ . All transitions are first-order. 
The surface field enhances the stability of the lower surface 
state, for which (j}{0) > and (j}{oo) < 0, and decreases the 
stability of the other. The leftmost segment of the horizontal 
line at h — represents the transition between the two bulk 
states; these states are not homogeneous in the presence of 
the surface field. The other segments describe bulk-driven in- 
stabilities of the surface states; for example, the lower surface 
state cannot exist for h > 0. 

surface first as simply terminating the bulk, and then in 
addition as supplying a surface field coupling linearly to 
the order parameter. 

The important new result of our analysis is the quite 
unexpected existence of solutions localized at the sur- 
face, solutions which exist even if the surface field is zero. 
These solutions are energetically favored for temperature 
and field values that are outside but not too far from 
the closed loop within which the modulated bulk state is 
stable. When the surface field differs from zero, the loop 
breaks apart (as shown in Figures 4 and 5). 

We now turn our attention to the applicability of our 
results to cholesteric liquid crystals in a field It is 
obviously desirable to estimate the conditions of chiral- 
ity, temperature, field, and surface interactions for which 
the surface states should be observable. To do this, we 
should examine the relationship between the variables of 
the theory and the experimental variables, by comparing 
the expressions for J-2 and in section II of this pa- 
per with the analogous expressions in section III-C of [|j . 
It is reasonable however, and far simpler, to expect the 
loop regions to scale by the same factors; this should be 
true independent of the strength of surface interactions. 
From Figures 1 and 2, the outer (surface-state) loop ex- 
tends over the range —4.5 ^ r < 1 while the inner loop 
extends over —1.2 ^ r < j, about a factor of four. Ac- 
cordingly, we estimate the surface-state region to be four 
times the size of the undulating-state region in temper- 
ature. From Q, the undulating state should occur for 
intrinsic pitches in the range of 1260 — 630 nm, at elec- 
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trie fields of the order of a few hundreds of kV/cm, or 
magnetic fields of roughly 40 T; the temperature width 
was estimated to be a few tenths of a degree. These 
conditions are very difficult to achieve and account for 
the fact that the undulating state has not yet been ob- 
served, although some groups plan to attempt the exper- 
iments. The surface-state region is expected to be one 
degree wide. Techniques sensitive to birefringence near 
the surface, such as Brewster-angle ellipsometry ]lO| , pl] |, 
may be able to detect the surface states. 

The surface states should appear in incommensurate 
systems where the modulated phase is driven by a neg- 
ative gradient-squared term. Conditions may be favor- 
able in magnetic Lifshitz-point materials like MnP, or in 
Langmuir monolayers or diblock copolymers. Other sys- 
tems in which a modulated phase is driven by a negative 
gradient-squared term are sodium nitrite and thiourea 
p2[ ; related systems are quartz and berlinite, but for 
these the modulated phase is two-dimensional. 
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